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ABSTRACT 



Aims. Galaxy cluster surveys based on the Sunyaev-Zeldovich effect (SZE) mapping are expected from ongoing experiments. Such 
surveys are anticipated to provide a significant amount of information relevant to cosmology from the number counts redshift distri- 
bution. We carry out an estimation of predicted SZE counts and their redshift distribution taking into account the current cosmological 
constraints and the X-ray cluster temperature distribution functions. Comparison between local and distant cluster temperature distri- 
bution functions provides evidence for an evolution in the abundance of X-ray clusters that is not consistent with the use of standard 
scaling relations of cluster properties in the framework of the current concordance model. The hypothesis of some evolution of the 
scaling law driven by non-gravitational processes is a natural solution to this problem. 

Methods. We perform a MCMC statistical study using COSMOMC, combining current CMB observations from WMAP, the SNIa 
Hubble diagram, the galaxy power spectrum data from SDSS and X-ray clusters temperature distributions to predict SZE cluster 
number counts. 

Results. Models reproducing well the X-ray cluster temperature distribution function evolution lead to a significantly lower SZE 
clusters number counts with a distinctive redshift distribution. Ongoing microwave SZE surveys will therefore shed new light on intr- 
acluster gas physics and greatly help to identify the role of possible non-gravitational physics in the history of the hot gas component 
of x-ray clusters. 
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1. Introduction 

Knowledge on the galaxy cluster population has greatly pro- 
gressed, thanks both to targeted observations and to systematic 
surveys carried by space X-ray facilities. The hot intra-cluster 
gas is also known to interact with the incoming photons of the 
Cosmic Microwave Background (CMB), leaving a specific fre- 
quency imprint known as the thermal Sunyaev-Zeldovich effect 
( |Sunyaev & Zeldovich 1972) l. This change of the sky brightness 
of the CMB can be written as a function of frequency and the 
!* ■ Compton parameter, y, proportional to the integrated gas pres- 
. £h ! sure a l° n g the line-of-sight of the cluster. For an individual clus- 
. ter, the integrated Compton parameter Y is the value of y inte- 
5_i ■ grated over the solid angle subtended by the clusterwhich angu- 
[ lar size happens to be nearly independent of redshift. This turns 
the SZE effect in a very effective probe of clusters at cosmolog- 
ical distances. Furthermore, SZE selection is very attractive in 
principle as the signal depends on an integral of the intracluster 
gas pressure, independently of its spatial distribution. 

The relevance of these SZE cluster surveys for cos- 
mological application has been outlined in the past (1994, 
Markevitc h et al. 19941 [Barbosa et al. 1996), motivating the on- 
set of key science projects for some of the most excit- 
ing experiments in the microwave regime. In particular, the 
Planck Surveyor ESA mission and the South Pole Telescope 
dRuhl et al. 20041 1 (SPT hereafter) are crucial experiences to ex- 
plore the SZE number counts and their cosmological signif- 
icance. With redshifts obtained from the optical and infrared 



follow-up, such samples have the potential to complement other 
cosmological probes and constrain the matter density parame- 
ter, £Im, the amplitude of matter fluctuations measured by the 
parameter cr^, the dark energy content of the Universe, and its 
equation of state w dLima & Hu 20071 1. Besides these important 
cosmological inferences, the SZE cluster cartography may also 
reveal the cluster formation processes out to high redshifts, an 
aspect which is further explored in the present letter. 



2. Cluster modelling 

The modelling of the clusters population and its evolution 
needs two key ingredients. The first one is an expression of 
the mass function and its dependence on cosmology. This was 
first theoretically attempted by Press and Schechter in 1974 
dPress & Schechter 1974| l. Since that time, numerical simula- 
tions have shown that the mass function follows a simple and 
nearly universal scaling relation with a dependence actually 
close to the initial Press and Schechter proposition, with devia- 
tions much below the precision level needed for present day ap- 
plications dTinker et al. 20 08). The second key ingredient is the 
relation between cluster mass and observable quantities. From 
scaling arguments (Kais er 19911 1. we can relate the temperature 
measured from X-ray to the mass: 



T x =A TM (hM) 2l3 \Cl 
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in which A(z, Qm) represents the density contrast (relative to 
the background density) by which clusters are defined, h is the 
Hubble parameter and Ajm describes the normalization of the 
Tx - M relation. Hydrodynamical simulations of cluster forma- 
tion have provided mass-temperature relations that actually fol- 
low the above relationship with dispersion of the order of 20%, 
while the normalization constant Ajm is uncertain depending on 
the gas physics ( |Kay et al. 2007) 1. The magnitude of the SZE is 
controlled by Y, the integrated Compton parameter depending 
on the gas mass and the average gas temperature: 



ACDM (£2 M =0.317, a 8 =0.842, /?=0.669) 



Y = KM g T g D~ 



(2) 



where D a is the angular distance, M g the total gas mass of the 
cluster, T g its (mass-)average temperature and K is a normal- 
ization depending only on physical constants. This relation, be- 
ing independent of the actual spatial distribution of the hot gas, 
makes the SZE signal an appealing proxy of the total mass. If 
the gas was isothermal in clusters, one would have T g = Tx- 
However, it is known that clusters are not isothermal with the 
temperature declining in outer parts (Vikh linin et al. 2005b fol- 
lowing a scaling law. We can therefore still assume that the gas 
temperature follows some scaling laws but with a different nor- 
mali zation T g = £,Tx- Using the observed gas fraction in clus- 
ters dSadat et al. 20051 I. a dark matter mass profile according to 
numerical simulation of cold dark matter (Navarro e t al. 199Tb 
and the observed temperature profile (Vikh linin et al. 20051 1 one 
can estimate £ ~ 0.6, a value which is quite uncertain, but in 
the following this number is left as a quantity to be determined 
from observations. With the above mass-temperature relation, 
the Y — M relation can be evaluated numerically: 



Y= 1.816.10^A m / B M 5/3 /! 8/3 Q 



a (z, a M ) 
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with fs as the baryonic gas fraction, D the dimensionless 
part of the angular distance and Y expressed in arcmin 2 . In this 
work, we take as the universal baryonic fraction corrected for 
some possible depletion T: fs = T^- . For a typical Coma- 
like cluster with about 10 15 solar masses, observed at z = 1, 
and using the standard scaling laws, the typical Compton distor- 
tion is Y is of the order of 2.3 x 10~ 4 arcmin 2 , corresponding at 
143GHz to a flux decrement with respect to the mean CMB flux 
of -21mJy, or a temperature fluctuation of -650/iK, well above 
both Planck and SPT sensitivity thresholds. Assuming the above 
relations, we can now compute the SZE source counts as well as 
their redshifts distribution, as expected in a Friedmann-Lemaitre 
cosmological model: theoretical number counts predictions of 
cosmic objects can be obtained from the mass function describ- 
ing the formation and evolution of objects over redshift. In this 
way, cluster counts dN(z) for clusters with intrinsic Y greater 
than a sensitivity threshold To are described as: 



dN(z) 
dz 



f +0 ° dV dn 

dM 

Jm (z) dz dM 



(4) 



where the mass threshold Mo(z) is obtained from equation 
0} dV is the volume element, accounting for the geome- 
try and expansion rate of the Universe, and jj-j is the clus- 
ter mass function, depending on the matter and energy den- 
sities of the universe and on the initial power spectrum of 
mass fluctuations. On the other hand, the very same SZE clus- 
ter population is a source of non-gaussian secondary CMB 
anisotropies (Cole & Kaiser 1988, Komatsu & Kitayama 1999 
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Fig. 1. Cluster temperature distribution function in the local and 
distant universe. Triangles are our estimation of the tempera- 
ture distribution function. The stars represent the estimated tem- 
perature distribution function at high redshift from the 400deg2 
sample restricted to the redshift range [0.4-0.6]. The continu- 
ous line corresponds to the predicted temperature distribution 
adjusted to fit CMB, SNIa, LSS and the local temperature dis- 
tribution function. The predictions for this model at redshift 0.5 
(the lower continuous line) are far in excess of the abundances 
estimated from observations. For visual comparison, the dashed 
lines are the corresponding predictions with an evolving mass- 
temperature relation, eq (6), with a — -1 (not fitted). 



IBond et al. 20051 ITaburet et al. 2010b . connecting the present 
large scale structure to the primordial CMB spectrum. This of- 
fers a way to pin down the very nature of the astrophysical evo- 
lution with redshift of the cluster Y — M relation. 



3. The X-ray picture 

However, since cluster masses are hardly measurable di- 
rectly, it is necessary to use a relation between mass with 
some other observables, including its dispersion, to derive 
constraints from the observed temperature distribution func- 
tion of X-ray clusters (Oukbir & Blancha rd 1992b . The con- 
straints read as a thin region in the parameter space {Ajm, 
erg) ( |Pierpaoli et al. 2003) , a known degeneracy between Ajm 
and cr 8 clearly identified and quantified. In the present analy- 
sis we use the standard A-Cold Dark Matter picture (ACDM) 
and Bayesian Monte Carlo Markov Chain (MCMC) estima- 
tions. We use the COSMOMC parameters estimation package 
with constraints provided by CMB data from WMAP 7 years 
(IJarosiket al. 2010b . Supernovae SN la data from the SDSS 
compilation sample ( Kessle r" ;t al. 2009b . matter power spec- 
trum P(k) estimation from SDSS LRG DR7 (IReid etal. 2010b 
and our estimation of cluster temperature distribution at differ- 
ent redshifts. For this, we use the local temperature distribution 
function of X-ray clusters from a local sample of clusters with 
fluxes above 2x 10 _1I erg.s _1 .cm -2 in the [0.1-2.4] keV band and 
with galactic latitude \b\ > 20 degrees obtained from the cluster 
X-ray data base BAX. This results in a sample comprising 48 
clusters. The theoretical temperature distribution function was 
derived from the mass function with the above defined mass tem- 
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perature relation, using the Sheth, Mo and Tormen mass function 
dSheth et al. 200 H . We can then constraint the six cosmologi- 
cal parameters of this vanilla model (Q«, Q.b, t, h, n, <x 8 ) and 
the normalization Ajm of the mass-temperature relation. As an 
illustration, we provide the theoretical temperature distribution 
function computed for the median model (the normalization be- 
ing Ajm = 7.31 keV, Figure[TJ, compared to the observed local 
abundance of X-ray clusters. Clearly, this model reproduces well 
the local abundance of clusters. A previous claim that the pre- 
dicted abundance of X-ray clusters in a standard concordance 
cosmological model exceeds the observed abundance at large 
redshift led us to compute the high redshift (z ~ 0.5) temperature 
distribution function of clusters in the best concordance model 
and compare it to the one inferred from the high-redshift sub- 
sample of the 400deg2 survey (Vikhl inin et al. 2 009). We find 
the standard scaling relation vastly overpredicts the derived X- 
ray cluster distribution function at high redshift. If we allow for 
some redshift evolution to compare the predicted abundance of 
X-ray clusters a mass temperature relation, we may rewrite the 
Tx - M relation as : 



T x =A TM (hM) 2/3 \n 



A(z, Q M ) 
' 178 



1/3 



(1 + z) 



Q+or) 



(5) 



Such non standard scaling law has been advocated 
dVauclair et al. 200"3l with a ~ — 1 in order to match the abun- 
dance of high redshift X-ray clusters in a concordance cos- 
mology in flux selected surveys. The predicted temperature 
distribution function at high redshift with non-standard evolu- 
tion of the mass-temperature evolution (using a power-law in- 
dex of a = -1 without further adjustment) is clearly pre- 
ferred over the non evolving one. Potential explanations for 
such evolution are likely to lie in non-gravitational heating of 
the intra-cluster gas, possibly tracing sources of cosmic stellar 
formation or/and additional complex physical processes within 
clusters, an area that is the subject of intense investigations 
dVoit 20051 IMcNamara et al. 20051 |Cattaneo & Teyssier 2007| 
Ida Silva et al."2009l 

4. The anticipated Sunyaev-Zel'dovich number 
counts 

An early attempt to estimate SZE number counts using the 
mass function and scaling laws as well as an estimation 
of their contribution to CMB fluctuations on small scales 
has been performed by Bartlet and Silk (1994). The use 
of SZE counts as a cosmological probe (Bar bosa et al. 1 996) 
was proposed around the early framework of Planck Surveyor 
Phase A studies, and soon followed by detailed investiga- 
tions dAghanim et al. 1997] Holder et al. 20001 |Kay et al. 2001 



Bens on et al. 2002"T Diego et al. 2002) further exploiting the po 



tential of SZE surveys counts as a cosmic probe. Such estima- 
tions were usually performed under the assumption of the stan- 
dard scaling law of the mass temperature relation. Since the SZE 
number counts are by definition the number of galaxy clusters 
brighter than a certain flux threshold S v , this reduces to the num- 
ber of galaxy clusters more massive than a mass threshold Mo. 
depending on the redshift z, from the Y - M relation. The conse- 
quence of the dispersion is dealt by assuming that in practice it 
can be approximated by a shift in the mass-temperature normal- 
ization (Blanc hard et al. 20001 1. In order to work out a consis- 
tent modelling of the clusters population and its evolution in the 
same picture, another MCMC analysis was performed including 



the redshift evolution power-law index a as an additional free 
parameter and adding the high redshift temperature distribution 
functions as a complementary observational constraint. For each 
of these models, we can compute the corresponding SZE source 
counts as well as their redshift distribution of the clusters, al- 
lowing a comparison of the importance of evolution effect to the 
statistical uncertainty inherent to the set of observations under 
consideration. Here, we evaluate expected high-sensitivity flux 
limits for both Planck and SPT as fiducial numbers. An impor- 
tant source of complexity, and possibly uncertainty in this topic 
is the accuracy of our knowledge of the experiment selection 
functions (Mel in et al. 20051 Uuin et al. 20071 1. We focus on the 
theoretical curves, leaving to the future works the full under- 
standing of the selection functions used to find clusters which 
will need a comprehensive modeling of the observational and 
detection processes, including the understanding of the known 
complexities of the missions algorithms to identify clusters and 
measure its fluxes. Here as fiducial criteria, we adopted three 
different sensitivity thresholds: Y = 2.10~ 3 (Case 1), corre- 
sponding to a secure lower limit for the sensitivity of Planck, 
Y = 4.10~ 4 (Case 2) for SPT or for an optimistic flux limit for 
Planck, Y = 4.10~ 5 (Case 3) which could be achieved if the se- 
lection function is well under control for experiments like SPT 
dRuhl et al. 20041) or AMIBA dHo et al 20091 . From the above 
discussion, at the brightest fluxes, source counts should be rel- 
atively free from resolution and/or selection effects, and allow 
some detection of the signature of evolution. Our final predicted 
counts and the associated uncertainties are shown in Figure 2: 
we give the predicted counts falling in the central 68% interval 
using cosmological models fitted to the local temperature dis- 
tribution function assuming standard scaling law for the mass- 
temperature relation; on the same figure, we provide expected 
counts for self consistent modelling of the cluster population, 
allowing evolution of the mass-temperature relation and includ- 
ing the high redshift temperature distribution function data as an 
additional constraint. It is clear that the integrated source counts 
show significant differences between the two cases. The possible 
lack of accurate knowledge of the Y - M normalization and of 
the selection function is likely to limit the interpretation of SZE 
integrated counts. However, the redshift distribution is defini- 
tively different, with a clear drop of the high redshift tail after 
z > 1 when the X-ray motivated evolution is taken into account. 
Because the SZE effect is sensitive to the actual average tem- 
perature of the gas, the observations in SZE will allow to clar- 
ify whether the kind of evolution needed for reproducing X-ray 
number counts affects the overall cluster gas or only the central 
region from which the X-ray emission comes. This is rather en- 
couraging as the high redshift clusters are likely to be point like 
in a moderately low resolution survey as will be provided by 
Planck, and thus alleviating the complex and possibly unsecured 
modelling of their spatial distribution for the determination of 
the selection function. 

Conclusion 

Current cosmological data and existing surveys of X-ray clus- 
ters provide tight constraints and thereby accurate predicted SZE 
number counts. The X-ray data probe an inconsistency between 
the theoretically expected and the observed high redshift tem- 
perature distribution function under the current vanilla ACDM 
cosmology if the standard scaling of the mass-temperature re- 
lation is used. The hypothesis of some evolution of the scaling 
law driven by non-gravitational processes is a natural solution 
to this problem. Our MCMC statistical study shows that evo- 



3 



P. Delsart, D. Barbosa & Alain Blanchard: SZ number counts 



Source counts 



Redshift distribution 



1000.0 



100.0 



"D 
T3 



CD 
CD 

CD 

■a 



S 10.0 




10 



CD 

§> 1 .000 



"O 

N 
T3 



n 0.100 



0.010 



10 



10' 



10" 



I I 

No evolution 
Evolution 

Sensitivity (tog to bottom) 
= 4\ 1 0* 
Yl=4A0t 

yf=2.io 



Yx£/0.57 (arcmin ) 




Fig. 2. Predicted galaxy clusters number counts detected from the Sunyaev-Zeldovich effect, after joint analysis of CMB, SNIa, 
LSS data and estimated temperature distribution function. Models describing self consistently the low and high redshift temperature 
distribution functions, needing non-standard evolution of Tx - M relation, lead to lower abundance of clusters, showing a deficit in 
the cluster number increasing with redshift (right). 



lution is required and that such evolution lead to a significant 
lowering of the anticipated SZE cluster number counts. Indeed, 
this may explain part of the statistical discrepancies suggested 
by South Pole dLueker et al. 2009b on the presence of a SZE 
signal anisotropy contribution to CMB fluctuations lower than 
expected in ACDM. This may also result from the geometri- 
cal structure of clusters (Ta buret et al. 20101 . While the X-ray 
data allowed the highlight of a redshift evolution in the galaxy 
cluster population, the SZE number counts have the potential to 
provide deeper insight on the actual nature of this evolution, by 
providing information allowing to test whether gas clusters did 
undergo a significant amount of non-gravitational heating affect- 
ing their global energy budget. This would call for further dedi- 
cated investigations in future analysis on cluster abundance and 
formation with potential consequences on the design of future 
experiments. 
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